Astronomy <fc Astrophysics manuscript no. 3187 


February 5, 2008 


(DOI: will be inserted by hand later) 





The limb darkening of a Cen B 
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Abstract. For the nearby dwarf star a Cen B (K1V), we present limb darkening predictions from a 3D hydro- 
dynamical radiative transfer model of its atmosphere. We first compare the results of this model to a standard 
Kurucz's atmosphere. Then we use both predictions to fit the new interferometric visibility measurements of 
a Cen B obtained with the VINCI instrument of the VLT Interferometer. Part of these new visibility measure- 
ments were obtained in the second lobe of the visibility function, that is sensitive to stellar limb darkening. The 
best agreement is found for the 3D atmosphere limb darkening model and a limb darkened angular diameter of 
#3D = 6.000±0.021 mas, corresponding to a linear radius of 0.863±0.003 Rq (assuming n = 747. 1± 1.2 mas). Our 
new linear radius is in good agreement with the asteroseismic value predicted by Thevenin et al. 1200211 . In view 
of future observations of this star with the VLTI/AMBER instrument, we also present limb darkening predictions 
in the J, H and K bands. 

Key words. Stars: individual: a Cen, Techniques: interferometric, Stars: binaries: visual, Stars: fundamental pa- 
rameters, Stars: atmospheres 



1. Introduction 

The limb darkening (hereafter LD) is a well known effect in 
stellar physics. Its manifestation is a non-uniform bright- 
ness of the disk whose edges appear fainter than the cen- 
ter. This effect occurs because of the decrease of the source 
function outwards in the atmosphere. The disk center then 
shows deeper and warmer layers whereas the edges show 
higher and cooler material. This means that the analy- 
sis of the intensity I\{n) at different latitudinal angles 
fi = cos 9 provides information on the temperature varia- 
tion with depth in the external layers of the star. This is 
therefore an excellent constraint to test atmospheric mod- 
els, validate or invalidate assumptions used to derive these 
models (like NLTE/LTE) and provide hints for improv- 
ing the input physics (equation-of-state and/or opacities 
in particular). The center-to-limb variation of the Sun is 
known for many years and has been measured for numer- 
ous ju and A (e.g. Pierce & Slaughter 1977, Neckel & Labs 
1994, Hestroffer & Magnan 1998) leading to a plethora of 
theoretical works which have improved our knowledge of 
the external layers of the Sun. 
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Traditionally, the analysis of solar and stellar LD 
is made by adopting an approximated law for 
generally a polynomial expansion in /i, either linear 
or non-linear (see e.g. Claret 120001 for recent develop- 
ments), whose coefficients are determined from ID atmo- 
spheric models, like ATLAS (Kurucz I1992f) or Phoenix 
(Hauschildt 1999). However, in spite of the detailed 
physics included in these codes, their ID nature is a lim- 
itation for deriving realistic emergent intensities. Indeed, 
these codes contain free parameters like the well-know 
mixing length parameter that are injected artificially in 
order to reproduce the properties of the turbulent convec- 
tion at the stellar surface. As a consequence, the compar- 
ison between these ID models and observations depend 
on the input parameters that thereby creates an impor- 
tant source of uncertainties. Moreover, the convection is 
by nature a 3D process. Its manifestation is the presence of 
bright granules and dark intergranular lanes. Reducing it 
to a ID process, i.e. ignoring horizontal flows and temper- 
ature inhomogeneities, changes the pressure scale height, 
the location of the surface and therefore also the emer- 
gent intensity (see e.g. Allende Prieto, Asplund & Fabiani 
Bendicho 2004|and Asplund et al. l20U0al for a comparison 
between multi-dimensional simulations). 
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The precise measurements of the center-to-limb vari- 
ation achieved nowadays require realistic stellar atmo- 
spheric models that take into account all the complexity of 
the stellar surface, and motivates the use of the new gen- 
eration of 3D radiative hydrodynamical (hereafter RHD) 
simulations. 

In this paper we propose a study of a Cen B 
(HD128621), a nearby K1V dwarf star. It is part of a visual 
triple star system whose brightest component, a Cen A 
(HD 128620), is a G2V dwarf. The motivation for the selec- 
tion of this star in the present work lies in both theoretical 
and observational considerations. From the interferomet- 
ric point of view, the proximity (1.3 pc) of the star is a rare 
opportunity to allow interferometric measurements since 
most of the nearby dwarfs have too small angular diame- 
ters to be measured. Our interest for this star has grown 
recently since our new measurements provide data points 
in the second lobe of the visibility function which is sensi- 
tive to the LD of the star. From a theoretical point of view, 
this star is important for various reasons. In particular, the 
recent detection of solar-like oscillations in a Cen A and B 
(Bouchy & CarrierUm I200T1 Carrier & Bourban I2005|l 
have led several authors (e.g. Morel et al. 120001 Thevenin 
et al. 1217021 Thoul et al. 12031 Eggenberger et al. 12004) to 
build evolution models of these two stars that are strongly 
constrained by the measured frequency spacings. The re- 
sult is a better, but still debated determination of the 
fundamental parameters of the system. 

In Sect. |21 we report the new interferometric measure- 
ments of a Cen B obtained since 2003 using the VINCI 
instrument. Sect. describes our 3D simulations to derive 
self-consistent stellar limb darkening of a Cen B. They are 
subsequently used to compute visibility curves in the near- 
infrared (SectQJ in order to interpret our measurements in 
terms of stellar angular diameter and to discuss the agree- 
ment between 3D limb darkening model and our second 
lobe visibility measurements. We also use our simulations 
to predict future observations (J, H, and K bands) that 
will be made with the next generation of instruments of 
the VLTI, such as the new AMBER instrument (Petrov 
et al. 120001 Robbe-Dubois et al. 12002)) . 

2. New interferometric observations 

A total of 37 new interferometric measurements of a Cen B 
have been obtained in 2003 on two baselines, D1-B3 (24 m 
in ground length) and B3-M0 (140 m), using the VINCI 
instrument (Kervella et al. 120001 Kervella et al. [2003a). 
The points obtained on the longer baseline are located in 
the second lobe of the visibility function of a Cen B, whose 
shape depends in the limb darkening. We have obtained 
1 000 interferograms on the B3-M0 baseline, in two series. 
Out of these, 534 were processed by the VINCI pipeline. 
The B3-M0 baseline observations are made difficult by the 
very low V 2 of the interferometric fringes, less than 2%. 
However, Fig. shows as an example the power spectral 
density of these very low visibility fringes where no bias is 
present. On the D1-B3 baseline, we have recorded 17500 
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Fig. 1. Average wavelets power spectral density (WPSD) 
of 299 interferograms of a Cen B obtained on JD 
2452770.6605 (11 May 2003). In spite of the very low vis- 
ibility (V 2 = 1.38%), the subtraction of the background 
noise (dotted line) from the processed fringes power peak 
(dashed line) leaves no residual bias on the final WPSD 
(solid line). The power integration is done between wave 
numbers 1970 and 7950 cm -1 . 

interferograms in 35 series (15141 processed). These new 
measurements were added to the V 2 values obtained on 
the E0-G1 baseline and already published in Kervella et 
al. I)2003bl hereafter Paper I) . The resulting squared visi- 
bilities are listed in Tables [2] and [3] 

We have used several stars from the Cohen et al. (1999 ) 
catalog as calibrators to estimate the point source re- 
sponse of the interferometer. They were observed im- 
mediately before or after aCenB. On the D1-B3 base- 
line, we have used HD 119193 (0 UD = 2.03 ± 0.022 mas), 
58Hya (6> UD = 3.13 ± 0.030 mas) and HD 112213 (0 UD = 
3.14 ± 0.025 mas). Approximately one third of the mea- 
surements was obtained with each of these calibrators. On 
the B3-M0 baseline, we have relied on HR4831 (#ud = 
1.66 ± 0.018mas), whose small size results in a low sys- 
tematic uncertainty on the calibrated V 2 values. The an- 
gular diameter estimates from Cohen et al. H1999[) have 
been verified by Borde et al. I|2002|l to be reliable within 
their stated error bars. The squared visibilities were de- 
rived using the processing methods described in Kervella 
et al. (|200 4i) . As an example, the calibration sequence used 
for the longest baseline B3-M0 is presented in Tabled 

3. Simulation of a 3D atmosphere 

In order to model the intensity profile of a Cen B, we have 
performed realistic, time dependent, three-dimensional ra- 
diative hydrodynamical simulations of its surface. The 
emerging intensity of the atmospheric model in different 
directions is used to built theoretical monochromatic limb 
darkening profiles for various wavelengths covering the 
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Table 1. Calibration sequence of aCenB on the B3-M0 baseline (140m ground length). The expected visibilities 
Vj-heo gi ven m this table include the bandwidth smearing effect. The interferometric efficiency given in italic characters 
corresponds to the value assumed for the calibration of these particular aCenB observations (see Tables |2 and [3J . 
The data of HR4831 marked with (*) were taken 2h before aCenB. They are listed to show the stability of the IE, 
but were not used for the IE estimation. 



JD 


Scans 


B(m) 


Azim. 


fi^± stat. (%) 


V$ MO ± syst. (%) 


IE ± stat. ± syst. (%) 


Target 


2452770.5474 


365 


139.309 


49.24 


29.37 ±0.39 


49.93 ± 0.80 


58.83 ±0.79 ±0.94* 


HR4831* 


2452770.5523 


316 


139.131 


50.51 


29.38 ±0.51 


50.06 ±0.80 


58.70 ± 1.02 ±0.93* 


HR4831* 


2452770.5572 


296 


138.913 


51.77 


30.05 ± 0.67 


50.19 ±0.80 


59.88 ± 1.34 ±0.95* 


HR4831* 


2452770.6368 


405 


128.934 


71.28 


32.05 ±0.35 


55.48 ± 0.75 


57.78 ±0.62 ±0.77 


HR4831 


2452770. 64f9 


408 


127.841 


72.49 


31.60 ±0.36 


56.01 ±0.74 


56.42 ±0.64 ±0.74 


HR4831 


2452770.6469 


392 


126.698 


73.69 


32.72 ± 0.36 


56.54 ± 0.73 


57.87 ±0.64 ±0.74 


HR4831 


2452770.6605 


299 


133.838 


59.85 


0.791 ±0.039 




57.36 ± 0.82 ± 0.74 


aCenB 


2452770.6656 


235 


133.277 


61.33 


0.777 ± 0.082 




57.36 ± 0.82 ± 0.74 


a Cen B 



Table 2. a Cen B squared visibilities. 
JD B (m) Azim. V' z ± stat. ± syst. (%) 



spectral domains of interest for the VINCI and AMBER 
instruments (in the 1.0-2.4^m range). 



Table 3. a Cen B squared visibilities (continued). 



JD B (m) Azim. V' 2 ± stat. ± syst. (%) 

E0-G1* 



2452462.5836 


60.441 


157.57 


17.02 


± 


0.36 


± 


0.26 


2452462.5870 


60.544 


158.40 


17.01 


± 


0.23 


± 


0.26 


2452462.5905 


60.645 


159.26 


16.80 


± 


0.77 


± 


0.26 


2452462.5949 


60.767 


160.35 


16.05 


± 


0.68 


± 


0.24 


2452465.6268 


61.541 


170.27 


16.76 


± 


1.05 


± 


0.26 


2452470.6203 


61.621 


172.05 


14.94 


± 


0.44 


± 


0.23 


2452470.6234 


61.650 


172.82 


15.59 


± 


0.42 


± 


0.24 


2452470.6278 


61.687 


173.92 


16.70 


± 


0.44 


± 


0.25 


B3-M0 
















2452770.6605 


133.838 


59.85 


1.379 


± 


0.07 


± 


0.02 


2452770.6656 


133.277 


61.33 


1.356 


± 


0.14 


± 


0.02 



* E0-G1 measurements reported by Kervella et al. J2003bt . 



3.1. The stellar atmosphere modelling 

The numerical code used for this work belongs to a new 
generation of 3D atmospheric codes developed for the 
study of solar (e.g. Stein & Nordlund ll1?gllM%jl and stel- 
lar (e.g. Nordlund & Dravins IT55U1 Asplund et al. IMJTH 
Allende-Prieto et al. l2~0"0"2l Ludwig et al. I30"0ljl g ranula- 
tion and line formations (e.g. Asplund et al. l2000al 12000 bl 
2000cll2IMll2nn5j) . The code solves the non-linear, com- 
pressible equations of mass, momentum and energy con- 
servation on an eulerian mesh. The 3D radiative trans- 
fer is solved at each time step along different inclined 
rays for which we have assumed local thcrmodynamical 
equilibrium (LTE). We have considered 10 latitudinal n 
points and 4 longitudinal ip points, and checked that a 
finer grid in (/x, ip) does not change the properties of the 
model. Realistic equation-of-state (including ionization, 
dissociation and recombination) and opacities (Uppsala 
opacity package, Gustafsson et al. I1975[) are used. The 
line-blanketing is taken into account through the opacity 
binning technic (Nordlund 1982). In the present simu- 
lation we have considered a cartesian grid of (x, y, z) — 
125 x 125 x 82 points. The geometrical sizes are 6 x 6 Mm 
for the horizontal directions and 5 Mm for the vertical one. 
The dimensions of this domain are large enough to include 
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a sufficiently large number of granules (n > 20) simultane- 
ously which prevents statistical bias. Periodic boundary 
condition is applied for the horizontal directions and trans- 
mitting vertical boundaries are used at the top and bot- 
tom of the domain. The base of the domain is adjusted to 
have a nearly adiabatic, isentropic and featureless convec- 
tive transport. The upper boundary is placed sufficiently 
high in the atmosphere so that it does not influence the 
property of the model. 

A detailed description of the current version of the code 
used in this paper may be found in Stein & Nordlund 
(1998). Unlike ID hydrostatic models that reduce all 
hydrodynamics into a single adjustable parameter, the 
present simulations are done ab initio by solving the com- 
plete set of RHD equations in a self-consistent way. All 
the dynamics and turbulence of the model come naturally 
from the equations of physics. Nothing is adjusted like 
the convective flux in the MLT formalism. The diagnostic 
made by such RHD simulations is therefore much more 
realistic than the ID models. We emphasize that the real- 
ism of these 3D simulations has been intensively checked 
for solar line formations (e.g. Asplund et al. 12000 bl l?000c 
2004), helioseismology (e.g. Rosenthal et al. 1999) and also 
for stellar line formations (e.g. Allcnde-Prieto et al. 2002). 

The adopted atmospheric parameters are those of 
Morel et al. pTOjl . i.e. T off = 5260 K, logg = 4.51 and 
[Fe/H] = +0.2. The simulation has been run for a few 
hours of stellar time which covers several convective turn- 
over times. The result is a three-dimensional, time depen- 
dent box representing the stellar surface. A snapshot of the 
disk-center surface intensity is represented on FigEl The 
structure of our model is similar to the one obtained by 
Nordlund & Dravins (1990) but is even more realistic since 
the present version of the code solves compressible equa- 
tions of hydrodynamics and uses more grid-points which 
allows a better treatment of the turbulence. 

3.2. 3D limb darkening 

The monochromatic surface intensity is computed for var- 
ious latitudinal [i and longitudinal if directions at the stel- 
lar surface. The limb darkening is obtained by hori- 
zontal (x, y), longitudinal and time averages of the surface 
intensity. For the time average we consider a sequence of 
2 hours of stellar time. The results are plotted in Fig. 
for the two extreme wavelengths of our spectral domain, 
1.0 and 2.4/im. For comparison, we overplot limb dark- 
ening obtained from a ID ATLAS9 model for the same 
wavelengths and for the same stellar fundamental param- 
eters. It appears that 3D RHD produces a less impor- 
tant center-to-limb variation than a ID static model. The 
departure from ID model increases with decreasing wave- 
lengths. Such behavior was also found by Allendc Prieto 
et al. (|2002|l for Procyon. However in the case of a Cen B, 
the departure from ID to 3D limb darkening is smaller, as 
a consequence of a less efficient convection in K dwarfs as 
compared to F stars. 




12 3 

X (Mm) 



Fig. 2. Snapshot of the disk-center (/i = 1) intensity 
emerging at the stellar surface at a representative time. 



The reason why the emergent intensity differs between 
ID and 3D models is due to the fact that the proper- 
ties of the superadiabatic and surface convective layers 
cannot be well described by the mixing length formalism, 
whatever parameter we choose. The temperature inhomo- 
geneities (granulation) together with the strong sensitivity 
of the opacity (H~) to the temperature make the warm 
ascending flows more opaque than they would be for an 
homogeneous ID model. This purely three-dimensional 
effect added with the contribution of the turbulent pres- 
sure, pushes the location of the surface to lower densities. 
The temperature gradient in these regions is steeper than 
in the ID case (see Nordlund & Dravins Il990|l . Since the 
continuum is formed in these layers, the emergent inten- 
sity is different. 

The correction due to 3D simulations (a few percents) 
is small but not negligible with respect to the precision 
obtained by the new generation of interferometric in- 
struments like VINCI or AMBER. This improvement is 
essential to derive an accurate angular diameter of the 
star. We report in Tabled our limb darkening predictions 
for a series of continuum wavelengths between 1.0 and 
2.4 /mi, corresponding to the JHK range accessible to the 
AMBER instrument. 



4. Visibility model and angular diameter of 
a Cen B 

In this Section, we describe the application of our 3D limb 
darkening models to the interpretation of the VINCI mea- 
surements of a Cen B. 
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0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

Fig. 3. Normalized monochromatic center-to-limb variation of the surface intensity obtained by 3D RHD 

simulations of a Cen B as a function of /i for two different wavelengths: 1 (im (left panel) and 2.4 /im (right panel), 
which correspond to the extreme limits of the spectral domain we have considered in this paper. The solid lines 
represent the 3D RHD limb darkening whereas the dashed lines represent limb darkening derived from ID ATLAS 
atmospheric models. In both case, the symbols o represent the values extracted from both ID (white) and 3D (black) 
simulations. 



Table 4. Limb darkening I(X, fi) for various wavelengths over the JHK range. 



A (H / n 


0.0 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


1.0 


1.050 


0.4434 


0.5745 


0.6453 


0.7069 


0.7605 


0.8087 


0.8527 


0.8932 


0.9311 


0.9667 


1.0000 


1.270 


0.4646 


0.6017 


0.6738 


0.7347 


0.7860 


0.8310 


0.8711 


0.9074 


0.9406 


0.9715 


1.0000 


1.650 


0.4838 


0.6752 


0.7487 


0.8039 


0.8462 


0.8812 


0.9110 


0.9369 


0.9601 


0.9811 


1.0000 


2.000 


0.5442 


0.7063 


0.7707 


0.8202 


0.8585 


0.8905 


0.9178 


0.9417 


0.9630 


0.9825 


1.0000 


2.200 


0.5729 


0.7220 


0.7817 


0.8283 


0.8646 


0.8950 


0.9211 


0.9439 


0.9645 


0.9831 


1.0000 


2.400 


0.5968 


0.7353 


0.7912 


0.8352 


0.8698 


0.8988 


0.9239 


0.9458 


0.9656 


0.9836 


1.0000 



4.1. Limb darkened disk visibility model 

In the simple case of a centro-symmetric star such as 
a Cen B, the visibility function measured using a broad- 
band interferometric instrument such as VINCI is defined 
by three wavelength dependent functions: 

1. The spectral energy distribution 5(A) of the star, ex- 
pressed in terms of photons (VINCI uses a photon 
counting detector). 

2. The wavelength dependent intensity profile of the star: 
7(A, M )//(A,1). 

3. The spectral transmission T(A) of the instrument, in- 
cluding the atmosphere, all optical elements and the 
detector quantum efficiency. 

Out of these three functions, T(X) is known from the con- 
ception of the instrument, as well as calibrations obtained 
on the sky (see Kervella et al. 2003b for details). The spec- 
tral energy distribution of the star S(X) can be measured 
directly using a spectrograph, or taken from atmospheric 
numerical models. 

From the 3D RHD simulations presented in Sect .[21 
we have obtained intensity profiles for ten distinct wave- 
lengths over the K band (chosen in the continuum). For 
each of these profiles, ten values of /i were computed. The 



resulting 10 x 10 elements 2D table /(A, /z)//(A, 1) was 
then interpolated to a larger 60 x 50 elements table, with 
a lOnm step in wavelength (over the 1.90-2.50 /im range), 
and a 0.02 step in fi. This interpolation preserves well the 
smooth shape of the intensity profile function. This proce- 
dure was also used to build the /(A, (i)/I(X, 1) table based 
on ID Kurucz model. The original sample (10 x 20) was 
interpolated to the same final grid as the 3D model. 

We can derive the monochromatic visibility law 
V\(B, 9) from the monochromatic intensity profile /(A, /j) 
using the Hankel integral: 

V x (B,9) = jj\(X,ti)J (^^Vl^^ 

where B is the baseline (in meters), 6 the limb darkened 
angular diameter (in radians), Jo the zeroth order of the 
Bessel function, A the wavelength (in meters), fi = cos 6 
the cosine of the azimuth of a surface element of the star, 
and A the normalization factor: 

A= f I(\,ii)ndii, (2) 
Jo 

To obtain the visibility function observed by VINCI in 
broadband, we have to integrate this function taking into 
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account the transmission of the instrument and the spec- 
tral energy distribution of the star: 



V K {B,B) = 



T(X) S(X)f A 2 dX 



J K [T(X) S(X)Y X 2 dX 



(3) 



Note the A 2 term that is necessary as the actual integra- 
tion of the squared visibility by VINCI over the K band 
is done in the Fourier conjugate space of the optical path 
difference (expressed in meters), and is therefore done in 
wavenumber a = 1/A. This corrective term ensures that 
the integration of the spectrum of the star is done precisely 
in the same way as in the instrument. 

This formulation is very general, as it does not make 
any particular assumption on the spectrum of the star, 
or on the wavelength dependence of its intensity profile 
7(A, M ). 

4.2. Fit of the interferometric data and angular 
diameter of a Cen B 

Considering the model discussed in Sect. 14. ll we now de- 
rive the limb darkened angular diameter #ld of aCenB. 
It is obtained by a standard x 2 analysis of the data. We 
define the reduced \ 2 of our fit by 



Xrcd 



N 



N 

E 



V 2 - V 2 A , 

% model 



(4) 



where n is the number of variables (n = 1 for our fit), 
TV the total number of measurements, i the index of a 
particular measurement, and cr; the standard deviation of 
the measurement with index i. 

The x 2 minimization was computed for three center- 
to-limb models: uniform disk (UD), ID ATLAS and 3D 
RHD. In each case the broadband square visibility curve 
Vk(B,9) is shown in Fig. 0] and Fig.0 In addition to 
the purely statistical error, we must also take into ac- 
count two systematic error sources: the calibration un- 
certainty and the wavelength uncertainty. The calibration 
uncertainty comes from the errors on the a priori angular 
sizes of the five calibrators that were used for the VINCI 
observations. It amounts to 0.012 mas on the final an- 
gular diameter. The wavelength uncertainty comes from 
the imperfect knowledge of the transmission of VINCI, in 
particular of its fluoride glass optical fibers. This trans- 
mission was calibrated on sky (Paper I), and the uncer- 
tainty on this measurement is estimated to be 0.15%. 
As it impacts linearly on the angular diameter value, it 
corresponds to 0.009 mas. These two systematic factors 
add up quadratically to the 0.013 mas statistical uncer- 
tainty, and result in a total error of 0.021 mas on the an- 
gular diameters of aCenB. The best fit angular diam- 
eter that we derive using our 3D limb darkening model 
is 03D = 6.000 ± 0.021 mas. The ID model produces a 
slightly larger diameter, dm = 6.017 ± 0.021 mas and 
the UD disk produces naturally a much smaller diame- 
ter, 6>ud = 5.881 ± 0.021 mas. 




50 100 150 

Projected Baseline (m) 



Fig. 4. Overview of a CenB squared visibilities. The con- 
tinuous line represents the broadband, limb darkened disk 
visibility model derived from the 3D RHD, with 6*3D = 
6.000 mas. 



There is no significant departure between the three mod- 
els in the first lobe of visibility. However, different ampli- 
tudes of the second lobe are observed. While UD model 
produces higher visibilities, the ID limb darkened model 
leads to slightly too low visibilities compared to our obser- 
vations. Overall, the 3D model leads to a slightly better 
agreement with observations. 

As expected, the difference 3D/1D is rather small since 
we are working in the near-infrared (K-band) and for a 
dwarf star. It is nonetheless comparable with o~ s tat, and 
therefore significant. 

4.3. Linear diameter 

Assuming the parallax value of Soderhjelm l|1999[) . tt = 
747.1 ±1.2 mas 1 , we find a linear radius of O.863±O.OO3R 
which agrees with results presented in Paper 1. We esti- 
mate that the adopted uncertainty in the a priori assumed 
(of T c fF = 50 K) leads to an error of about 0.0003 Rq, i.e. 
ten times smaller than the derived uncertainties. From the 
ID analysis, we derive a radius of 0.865 ± 0.003 R Q , larger 
than the radius found by the RHD approach by about 
lostat. In addition to the corrections it provides, the use 
of 3D simulations is also motivated by the absence of ad- 
justable parameters, which is not the case for ID models. 

This slightly smaller linear radius obtained from 3D 
RHD simulations compared with the one derived from ID 
ATLAS model supports the suggestion of a smaller mass 



1 One should note that there is a rather broad distribution 
of parallax values for a Cen in the literature, as discussed in 
Paper I. While the value from Soderhjelm ! 1999: is the result 
of a careful reprocessing of the Hipparcos data, the possibil- 
ity of a bias beyond the stated 1 a error cannot be completely 
excluded, in particular due to the extreme brightness and bi- 
narity of a Cen. 
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Fig. 5. Close up views of the squared visibilities of aCenB in the lower part of the first lobe (left panel) and the 
second lobe (right panel). The continuous line represents the broadband, limb darkened disk visibility model derived 
from the 3D RHD, with 9^ = 6.000 mas. The dashed lines correspond to results obtained from ID ATLAS model 
with 6*id = 6.017 mas. The upper dotted curve is a UD model with #tjd = 5.881 mas. 



(M = 0.907 M Q , Kervella et al. 2003) rather than the 
larger one (M = 0.934 ± 0.007 A/©) proposed by Pourbaix 
et al. I|2002|l . However, stellar evolution models are sensi- 
tive to many parameters and a smaller radius do not al- 
ways lead to a smaller mass. More investigations are thus 
needed before we can reach a definite conclusion about 
the mass of a Cen B. In this context, our improved radius 
provides an additional constraint on asteroseismic diag- 
nostics. 

5. Conclusion 

In this paper we improve the determination of the radius of 
a Cen B in two respects. Firstly, we report the first inter- 
ferometric measurements in the second lobe of visibility. 
Secondly, in order to derive a reliable value of the angu- 
lar diameter of the star, we have performed realistic 3D 
RHD simulations of the surface of a Cen B. By compari- 
son with observations we find a radius of 0.863 ±0.003 Rq. 
The correction provided by the 3D approach is small but 
significant (especially in the K band probed with VINCI) 
since it provides a radius smaller by roughly 1 a s t a t com- 
pared with what can be obtained by ID models. Moreover, 
the use of 3D RHD is preferable since it does not intro- 
duce adjustable parameters to describe convection. We 
also emphasize that for hotter A-F stars the correction due 
to 3D analysis will be larger than for a Cen B. Though it 
is small, we have shown that even for a K-dwarf like a Cen 
B, the correction obtained by the use of RHD simulations 
should not be neglected for very high precision interfero- 
metric measurements. In the next few years, the combi- 
nation of high visibility precision and long baselines will 
require the use of realistic theoretical models of the stellar 
limb darkening to extract the true photospheric angular 
diameter of the observed stars from the observed visibil- 
ities. Conversely, observations beyond the first minimum 



of the visibility function will sample directly the light dis- 
tribution on the surface of the stars, therefore providing 
constraints for the atmosphere structure models. Future 
observations with the VLTI will allow to sample tightly 
the second lobe of the visibility function of several solar 
analogs (including a Cen A and B), and therefore derive 
their intensity profiles with high accuracy. The compar- 
isons between our theoretical predictions of limb darken- 
ing and the future observations made by AMBER will 
be an excellent test for our modelling of the surface of 
a Cen B. Indeed, AMBER will provide new interferomet- 
ric data which will contain much more information com- 
pared with VINCI. There will be two major advantages 
with AMBER: 

— It will provide a wavelength dependence of the visibil- 
ity ([1.9-2.4] /im) therefore allowing differential studies 
of the limb darkening as a function of wavelength. 

— AMBER can combine simultaneously the light from 
three telescopes and therefore measures the closure 
phase. This gives an advantage to determine the an- 
gular size of the star when observing in the minima of 
the visibility function. 

These improvements will lead to better constrained an- 
gular diameters of a Cen A and B, and therefore a high 
precision measurement of the ratio of the linear radii of A 
and B, independent of the parallax. 
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